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ABSTRACT 

Context. Providing a theoretical basis for the baryonic Tully-Fisher Relation (BTFR; baryonic mass vs rotational velocity in spiral 
galaxies) in the ACDM paradigm has proved problematic. Simple calculations suggest too low a slope and too high a scatter, and recent 
semi-analytic models and numerical galaxy simulations typically fail to reproduce one or more aspects of the relation. Furthermore, 
the assumptions underlying one model are often inconsistent with those behind another. 

Aims. This paper aims to develop a rigorous prediction for the BTFR in the context of ACDM, using only a priori expected effects 
and relations, a minimum of theoretical assumptions, and no free parameters. The robustness of the relation to changes in several key 
galactic parameters will be explored. 

Methods. I adopt a modular approach, taking each of the stand alone galaxy relations necessary for constructing the BTFR from 
up-to-date numerical simulations of dark halos. These relations - and their expected scatter - are used to describe model spirals with 
a range of masses, resulting in a band in the space of the BTFR that represents the current best guess for the ACDM prediction. 
Results. Consistent treatment of expected ACDM effects goes a large way towards reconciling the "naive" slope-3 ACDM prediction 
with the data, especially in the range 10 9 M o < M\, M < 10" M Q . The theoretical BTFR becomes significantly curved at Mt, ar > 10" M Q 
but this is difficult to test observationally due to the scarcity of extremely high mass spirals. Low mass gas-rich galaxies appear to 
have systematically lower rotational velocity than the ACDM prediction, although the relation used to describe baryon mass fractions 
must be extrapolated in this regime. The fact that the BTFR slope derived here is significantly greater than in early predictions is a 
direct consequence of a corresponding increase in the expected sensitivity of baryon mass fraction to total halo mass. 

Key words. Galaxies: fundamental parameters - Galaxies: spiral - Galaxies: halos - dark matter 
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1.1. The Tully-Fisher Relation in ACDM 

The Tully-Fisher Relation (TFR) was originally proposed as a 
correlation between 21 cm line width and optical luminosity in 
spiral galaxies ( Tully & Fisher |1977| l, but it has subsequently 
become accepted that these are observational proxies for more 
fundamental physical properties, namely rotational velocity and 



stellar or baryonic mass (e.g. Verheijen 2001 McGaugh et al. 



2000). The canonical approach towards deriving the TFR adopts 
the ACDM model of cosmology, in which the universe is flat and 
consists (by mass) of 4.6% baryons, 23% cold dark matter, and 
73% dark energy ( Komatsu et al.|2009) . In this paradigm, visi- 
ble galaxies are surrounded by roughly spherical halos of dark 
matter. Since it is total enclosed mass (baryonic plus dark) that 
determines a galaxy's rotational velocity at any particular radius, 
it is clear that the TFR describes a coupling between the baryons 
in a galaxy and the surrounding dark matter. The extent to which 
such a coupling is expected for cold, weakly interacting dark 
matter has recently been a hotly-debated subject, and consensus 
has not yet been reached (e.g. Foreman & Scott|2012; McGaugh 



|201 la) . The situation is complicated by the fact that many differ- 
ent forms for the TFR exist in the literature, using luminosities 
(at various wavelengths), stellar mass or baryonic mass as the 
dependent variable, and one of several different rotational veloc- 



ity measures as the independent variable (e.g. see Avila-Reese 
|et al.|2008[ [Masters et al.|2008) . 

The earliest ACDM predictions for the TFR started with the 
expected relation between the total mass of a galaxy and its 
halo, Myir, and the corresponding characteristic velocity disper- 
sion V v h , assuming virial equilbrium up to a radius R vn . Within 
this, the galaxy's average mass is taken to be a fixed multiple of 
the background mass density of the univ erse; simple manipula- 
tions then yield M vil oc ( White|l997|l. However, M vil and V vir 
cannot be measured directly: tfiey are theoretical quantities per- 
taining to the putative dark matter halo. In order to make contact 
with the observed TFR, one must convert halo mass to observed 
mass or luminosity, and virial velocity to disk rotational velocity. 
The expectation was that the baryonic mass of a galaxy would be 
a fixed fraction of its total mass ( |White||1997[ |Mo et al.|1998) , 
equal to the universal ratio of baryonic to total mass (now be- 
lieved to be 0.17; Komatsu et al. 2009). Assuming further that 
galaxies' rotational velocities are proportional to their virial ve- 
locities, one derives a prediction for the baryonic mass TFR 
(BTFR): Mb a r °c V^ ot . Such reasoning would provide a conve- 
nient explanation for a BTFR with a slope of 3 in log-log space. 
However, observational studies persistently find a slope greater 
than 3, typically in the range 3.5-4 (e.g. see Torres-Flores et al 
|20TT] table 3). 

Several effects might be expected to invalidate the premises 
used to derive the slope-3 BTFR: 
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There is no a priori reason why a galaxy's rotational veloc- 
ity should be proportional to its virial velocity. A better ap- 
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proximation would be to use the density profile of a galaxy 
and halo to compute the overall rotation curve, and to take 
the velocity at some characteristic point on this as the TFR's 
independent variable. The maximum rotation velocity will 
be used in this paperQ For example, dark halos are com- 
monly parametrized by the NFW density profile (Navarro 
|et al.|1996> : 



P(r) = Pcrit 



ZL(1 + _L)2' 



(1) 



(Mo 



with r s the scale-length and 60 a characteristic overden sity. 
<So may be written in terms of the concentration, c ■■ 
|et al.|1998> : 

A r 3 
So= - 



3 ln(l + c) - 



(2) 



l+c 



A is the mean density of the galaxy, at the time of virialisa- 
tion, in multiples of the critical density of the universe, p cr i t . 
It will be set to 200 in this paperj^J 

The ratio of maximum rotational velocity to virial velocity 
depends on the concentration (Bullock et al. 2001), which 
is found in N-body simulations to be a function of virial 
mass ( jGao et al.|[2008| |Zhao et al.|[2009] ). After changing 
independent variable from virial velocity to maximum ro- 
tation velocity, this M v j r -concentration relation is therefore 
expected to change the TFR's slope. 

2. So far, no mention has been made of the collapse of the 
baryons in galaxies into disks, yet both variables in the 
BTFR are properties of these disks. In ACDM, the parti- 
cles in a newly-virialised galaxy acquire angular momentum 
through the action of cosmological torques (Peebles 1969). 
Dissipative processes cause the baryons to lose energy, and 
hence fall towards the centre, and their rotation forces them 
into a flattened disk. This modifies the rotation curve of the 
galaxy. Furthermore, the distribution of dissipationless dark 
matter is affected by baryonic collapse in a process known as 
adiabatic contraction (AC). The shells that make up the dark 
halo are pulled gravitationally towards the disk, with their 
final positions fixed by conservation of specific angular mo- 
mentum dBlumenthal et al.|1986||Gnedin et al.|2004||20TT) . 
This effect is also expected to modify the galaxy's rotation 
curve, and hence the velocity used in the TFR. 

3. The proportion of a halo's mass that is in baryons is un- 
likely to be universal. In particular, baryons are expected 
to be more easily ejected (e.g. due to supemovae or stellar 
winds) from lower mass galaxies with shallower potential 
wells (e.g. |Gnedin|2TX)0"t|Hoeft et al.|2006| ). This would make 
the baryon mass fraction a rising function of total (virial) 
mass. As well as directly determining the baryonic mass that 
constitutes the dependent variable of the BTFR, this M v ; r - 
Mb ar relation governs the relative contributions of the disk 
and halo to the overall rotation curve, and hence affects the 
position of its maximum. 

Whilst the first two of these effects must be studied through 
numerical simulation of dark halos, the third is accessible ob- 
servationally via weak lensing and satellite kinematics (e.g. 

1 Selection effects arising from the use of the flat part of the rotation 
curve in some observational studies will be explored in Sect. 



4.3 



2 A determines the notional extent of the dark halo (R VK ), yet no con- 
sensus exists regarding exactly what value it should take (e.g. see |Zhao| 



|et al.|[2009 p. 2). Halos do not have a well-defined edge, but rather 



Mandelba um et al.|2006|[Conroy et al.|2007| ). These techniques 
do indeed suggest that M^/M^ rises with M v ; r , but may only 



be used for galaxies with M stal > 10 M Q (e.g. see Behroozi 
|et al.|[20T0] fig. 11) and must therefore be extended to lower 
mass to make contact with the majority of TFR measurements. 
A recent technique for doing this is called Halo Abundance 



Matching (HA M; Rodrfguez-Puebla e t al.|201lj |Behroozi et al 



|2"0T0l|Moster et al.|2010||Guo et al.|2010| ). The premise of HAM 



is that the most massive halos harbour the most luminous galax- 
ies, and thus galaxies observed in a particular survey may be 
uniquely assigned to dark halos formed in N-body simulations. 
In association with appropriate mass-to-light ratios for convert- 
ing luminosity to stellar mass, and a relation between gas mass 
and stellar mass, the output of HAM is an M v j r -Mb ar relation for 
the galaxies in the survey. Encouragingly, this technique agrees 
well with the direct methods at high M stal (e.g. see Behroozi et al. 
20101 fig. 11). 

The way in which AC and the M^-M^ and M v j r - 
concentration relations modify the "naive" BTFR prediction is 
the primary concern of this paper. Before discussing the gen- 
eral approaches that have previously been employed in studies 
of this kind, it is important to describe an additional tension be- 
tween observation and ACDM prediction that is exacerbated by 
the three effects listed above. It has been noted that the observed 
BTFR is very tight (e.g. Verheijen"]|2001| |Trachternach et a! 
2009 Mc Gaugh|[20T2) ). Indeed, McGaugh argues in |McGaugh 



( |2011b[ ) that the scatter of individual spirals from a power- law 
TFR may be entirely accounted for by observational uncertain- 
ties, suggesting that the theoretical relation has zero intrinsic 
scatter. However, several factors are expected to create signif- 
icant scatter in the ACDM prediction, including variations in 
the density profiles of dark halos caused by different mass ag- 
gregation histories (manifest in scatter in the concentration of 
halos of a given mass; e.g. see Eisenstein & Loeb]|1996| Jing 



2000) and in baryon mass fractions (e.g Behro ozi et al.| |2010; 
|Mo et alf| l998 ). These sources of scatter are described in more 
detail in Sect. 4.2 McGaugh argues that such effects will in- 



evitably make the predicted scatter irreparably larger than that 
in the data, leading him to the conclusion that the BTFR is 
evidence against the ACDM paradigm itself. McGaugh pro- 
poses that the BTFR data are more consistent with the expec- 
tations of Modified Newtonian Dynamics (MO ND), which im 



plies a BTFR of the form M bar oc V r 4 ot (pVlilgroml 1 983| iMcGaugh 



2011b). This appears to be in agreement with recent observa 
tions of gas-rich spiral galaxies in addition to many observations 
of star-dominated spirals (McGaug h|2012 hereafter MG12). 



1.2. The paths to an improved prediction 

Two main approaches exist for improving the "naive" TFR pre- 
diction by taking into account the effects described above. The 
first is to perform a full numerical simulation of one or more 



spiral galaxies (e.g. I Piontek & Steinmetz|2009| |Governato et al. 
2007 1 Tissera et al.|2010| ). Cosmological parameters and initial 



merge continuously with the background mass density of the universe. 



conditions are specified, and galaxies evolved stepwise through 
time according to general relativistic equations of motion and 
prescriptions for gas cooling, star formation, and stellar feed- 
back. The baryonic or stellar masses of the resulting galaxies 
can then be measured in addition to their rotational velocities, 
enabling their positions on the TFR to be determined. In the fu- 
ture, this approach will likely provide the most complete and 
robust estimate for the TFR expected in ACDM. 

Currently, however, such simulations face several problems 
that limit their usefulness, including insufficient resolution for 
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successful modelling of galaxy structure, insufficient comput- 
ing power for the simulation of a statistically significant num- 
ber of spirals, and uncertainties in our theoretical understand- 
ing of dissipative baryonic physics. These problems are mani- 
fest in significant discrepencies between simulated spirals and 



both the observed TF R (|Dutton et al 20 10| l and expected galaxy 
formation efficiencies ( |Guo et al.|2010) >. Furthermore, this holis- 
tic approach risks masking the effects that individual correc- 
tion factors have on the theoretical TFR: cosmological param- 
eters are essentially fed into a 'black box' describing a whole 
gamut of physical processes known to varying degrees of ac- 
curacy. Finally, approximations and uncertainties are introduced 
into the simulation that are not strictly required for prediction of 
(some versions of) the theoretical TFR. For example, star forma- 
tion rates and thresholds are necessary for determining the stellar 
mass and luminosity of a galaxy - and hence its position on the 
luminosity or stellar mass TFR - but not the total baryonic mass. 
As another example, consider the M v ; r -Mb a r relation. This may 
be predicted in numerical simulations from a prescription for 
the extent of baryon expulsion from a galaxy by supernovae and 
stellar winds, but may also be measured observationally (as doc- 
umented in Sect. [TTTi. A disagreement would suggest a flaw in 
the simulations which would propagate into the predicted TFR, 
but the specific effect of this would not be readily visible. Using 
instead the observed M v i r -M\, ar relation (or an extension thereof, 
such as that from HAM) would provide a more accurate pre- 
diction for the TFR, and would disentangle the result from any 
potential conflict between observed and simulated baryon mass 
fractions. 

The alternative to complete cosmological simulation is semi- 
analytic modelling (see e.g. Baugh 2006). In this modular ap- 
proach, the relationships between variables important for the 
TFR are determined empirically, where possible, and are oth- 
erwise given simple forms in accordance with the results of N- 
body simulations of dark halos. This allows the development of 
a TFR prediction with a minimum of assumptions and enables 
the individual effects of different components of the prediction 
to be investigated, thus helping to isolate the specific factors re- 
sponsible for potential disagreement with observation. The dis- 
advantage is that the various parameters used (those describing 
the Myir-Afbar and -concentration relations, adiabatic con- 
traction, and the halo density profile) are typically derived from 
different numerical simulations, which may use different initial 
conditions or halo virial parameters and therefore be marginally 
inconsistent. One of the first analyses of this nature was per- 
formed by |Mo et aT71 < |1998| hereafter M98), whose methodolog- 
ical framework provides a template for connecting the various 
pieces required for a prediction of the TFR. I adopt this template 
here. 

When the ingredient relations are not believed to be well- 
known, free parameters are often introduced which may be tuned 
to produce a desired result. An example of such an approach 



is|Dutton et al.|(2007 hereafter D07), where agreement with an 



observational TFR data set is optimised by searching through a 
high-dimensional parameter space for the region that minimises 
a x 2 goodness of fit estimator. This method is useful for telling 
us what the existence of the TFR implies for the properties of 
spiral galaxies. However, one cannot be said to have explained 
the TFR if one has simply selected parameter values that repro- 
duce it. A basic tenet of my methodology (which is described in 
more detail in Sects. [T~4| and |2|i is that the introduction of free 
parameters in the context of the TFR is not necessary and can in 
some circumstances produce misleading conclusions. The rea- 



sons for this, drawing examples from the methods and results of 
D07, are as follows: 



It is necessary to assume some simple form for the relations 
in the first place in order to specify them with only a few pa- 
rameters. However, these parametrizations themselves may 
lack physical motivation. For example, D07 use a power-law 
to describe the relation between and M v „ (a common 
approximation since the work of |Mo et al.||1998"] l, but re- 
cent Halo Abundance Matching studies observe a turnover 
at high masses ( Rodrfgu ez-Puebla et al.|20"TT Moster et al. 
2010). These results will be inconsistent with those of D07 
whatever parameter values are used. 

D07 tune their parameters to produce agreement with a 
power-law that is taken to fully describe the TFR. Yet no 
reason for adopting a power-law parametrization is inherent 
in the data - logically, such a form must instead be a prop- 
erty of a theoretical prediction. Indeed, Trujillo-Gomez et al. 
(201 1 hereafter TGI 1) find agreement between the data and 



a theoretical line that is significantly curved, suggesting that 
it is not necessary for a prediction to be a power-law in order 
to be acceptable. Furthermore, the slope and intercept of the 
power-law fit to the TFR may not be very well defined: it 
may be possible to fit the TFR almost equally well using two 
power-laws with totally different slopes and intercepts (e.g. 
Foreman & Scott|2012| fig. 2). An observed TFR comprising 
data points with uncertainties is not uniquely equivalent to a 
power-law with uncertainties in its slope and intercept. 
Even if a power-law could fully describe the observed TFR, 
its parameters are obviously dependent on the data set used. 
There have been many studies into the TFR, all yielding 
slightly different results. Thus the results in D07 are strictly 
valid only for the particular observational data that they use 
to construct the TFR, intertwining theory and observation in 
a non-trivial way. In fact, different best fit slopes and inter- 
cepts can be given even for the same data set, by the use of 
different fitting methods (e.g. see MG12, Sect. 2.6.1). Taking 
a particular power-law to be the correct description of the 
TFR requires the assumption that a particular fitting method 
is superior to all others. This is not the case: they are all 
just different approaches to an underdetermined statistical 
problem. It is not clear how changes to the parameters of 
the power-law TFR would affect the optimum values of the 
tunable parameters. 

Although it will be possible to find the global minimum of 
the x 1 value of the TFR fit, it may be that distant local min- 
ima exist in the many-dimensional parameter space with al- 
most as low a x 2 value. Given the uncertain nature of the 
observed TFR, these could provide equally good descrip- 
tions of disk galaxies and may, in fact, be better physically- 
motivated. 

All of the free parameters are correlated in the context of the 
X 2 value that they produce, making it difficult to assign them 
unique and independent uncertainties. 
The relationships obtained after tuning the free parameters 
lack physical significance because they have become di- 
vorced from the motivations that originally existed for intro- 
ducing them. At the end, the results must be compared with 
independent measurements to check consistency, but, to the 
extent that such independent measurements exist, one may 
as well use them from the outset. 
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1.3. Two important methodological issues 

In this section I discuss two further methodological issues that 
influence the nature of the TFR and the way in which it is pre- 
sented. By clarifying them here I intend to illuminate the ap- 
proach to TFR prediction that I believe to be most constructive. 



1 .3.1 . Different versions of the Tully-Fisher Relation 

As alluded to in Sect. |1.1| several different versions of the TFR 
exist. In particular, one may use luminosity, stellar mass or bary- 
onic mass as the dependent variable, and the significance that 
one attributes to these different TFRs depends on the perspec- 
tive brought to the problem. McGaugh, for example, believes 
the baryonic TFR to be fundamental because of its power-law 
nature and low intrinsic scatter - and because it is such a rela- 
tion that arises naturally in MO ND (|McGaugh||201 lb| MG12). 
Others (for example TG11 and |Dutton et al.|201 0) consider the 
TFR to be a consequence of complicated physical processes oc- 
curing during hierarchical galaxy formation. As such, the TFR 
is no more significant than the luminosity TFR (LTFR). Since 
it is luminosity that is directly observed, these authors take the 
latter as the touchstone for their ACDM models. (The baryonic 
TFR is clearly of secondary importance in TG11, where its in- 
trinsic scatter is not even considered.) However, regardless of the 
MOND-ACDM debate, it seems intuitively reasonable that mass 
should be more fundamental than luminosity when it comes to a 
relation with rotational velocity; luminosity depends on the de- 
tails of stellar populations, which are not relevant for determin- 
ing how fast a galaxy spins. Even in ACDM, one would expect 
M v j r to be correlated more strongly with M\, ar than M star , and 
hence the BTFR is again the more "fundamental". 

If one starts with a Myir-^bar relation yet wishes to construct 
a theoretical LTFR (as in D07, for example), two additional in- 
gredients are required. First, a prescription for star formation is 
neccesary to convert baryonic mass to stellar mass. This is not re- 
quired for predicting the BTFR - where stellar (and gas) masses 
are measured directly to constitute the observational data points 
- and unnecessarily introduces an extra degree of uncertainty. 
Put another way, comparison of the theoretical and observed 
LTFRs tests the hypothesis ACDM + star formation prescrip- 
tion, whilst comparing BTFRs tests ACDM alone. Disagreement 
with the expected LTFR may simply be due to errors in the star 
formation rate. The second required ingredient is a mass-to-light 
ratio to convert stellar mass to luminosity. Such a ratio is also 
required for consideration of the BTFR, but here for the obser- 
vations as opposed to the theory. Reducing the number of as- 
sumptions and uncertain parametrizations that go into the theory 
allows it to be more directly compared with a range of observa- 
tional data sets using different mass-to-light ratios. 

Nevertheless, there do exist uncertainties inherent in the 
BTFR that are absent from the LTFR. If gas masses in spirals 
are not well known, the BTFR may be subject to large system- 
atic uncertainties that render impossible a fair comparison with 
any ACDM prediction. In this vein, [Gnedin (2012) has argued 
that disagreement between BTFR observations and the ACDM 
prediction implies the existence of entirely unobserved baryons 
in spirals, which he believes to be in the form of warm ionised 
gas created by radiation from the inner stellar disk and the ionis- 
ing cosmic background. Such claims are diametrically opposed 
to those of McGaugh, who argues that systematic errors such as 
missing gas mass should not be able to produce a tight BTFR 
correlation where none would otherwise exist. In other words, 
McGaugh uses his zero-intrinsic-scatter slope-4 BTFR (MG12, 



fig. 1) as evidence that all baryonic mass has been accounted 
for: it would be extremely unlikely for omission of one baryonic 
component to turn the messy relation predicted by ACDM into 
the precise power-law relation that is observed. McGaugh also 
cites observational studies which find insufficient ionised gas to 
significantly affect the BTFR (e.g. Anderson & Bregman|2010) . 
In fact, it is argued in McGaugh ( 2011b| l that the uncertainties 
involved in determining the gas mass (due to imperfectly-known 
distances and uncertainties in measured line fluxes) are dwarfed 
by those associated with the conversion of luminosity to stel- 
lar mass (which stem from stellar population modelling and the 
Initial Mass Function). This suggests that of all the variants of 
Tully-Fisher Relation, the BTFR should provide the cleanest test 
of any theory. In any case, it is the BTFR that we must investigate 
if we are to make contact with McGaugh' s claim that the slope 
and scatter of this relation constitute evidence against ACDM. 

One final point has great significance in the context of my 
work. Even if much ionised gas were to exist, this would affect 
neither the BTFR data points nor a theoretical prediction using 
baryon mass fractions obtained by Halo Abundance Matching. 
This is because the baryonic mass used in HAM is precisely the 
baryonic mass that we observe, and thus identical to that plotted 
in the BTFR. Use of an empirical M v ; r -Mb al - relation from HAM 
renders the issue of BTFR prediction entirely independent of the 
missing mass problem (observed baryon mass fractions signifi- 
cantly different from the "average" cosmological value of 0.17). 
The conjecture of Gnedin ( 2012| l may shed light on the where- 
abouts of this missing mass, but has no ramifications for a BTFR 
built using HAM. 



1.3.2. Inclusion of elliptical galaxies 



Both TG11 and Dutton et al. ( |2010 1 introduce ellipticals into 



their TFR, thereby turning it into a more general "luminosity- 
velocity" relation. By enlarging the scope of the "TFR" in this 
way, they are able to check their predictions against more ob- 
servational data and explore the differences in the properties of 
early- and late-type galaxies. However, it is not clear to what ex- 
tent the direct comparison of spirals and ellipticals in this way 
is appropriate. V3cr/ os is used as an elliptical equivalent of ro- 
tational velocity, where cr 2 los is the stellar velocity dispersion 
along the line of sight to the galaxy. This assumes virial equilib- 
rium and an isotropic stellar velocity distribution; the (poorly- 
known) geometry of ellipticals can cause differences of up to 
20% ( |McGaugh & Wolf|[20T0) l. Further, one must be careful 
to measure <x/ os at a radius comparable to the one at which the 
rotational velocities of spiral galaxies are measured. The intro- 
duction of elliptical galaxies clearly creates several uncertainties 
that need not be of concern if one is dealing specifically with the 
TFR. 

TG11 observe a difference in the position of spiral and el- 
liptical galaxies in their mass-velocity plot (their fig. 11). As 
they themselves note, this is likely to be due to a dependence 
on galaxy morphology of the relationships which act as ingredi- 
ents of the theoretical BTFR, for example the M v ; r -Mt, al relation. 
No distinction is made between spirals and ellipticals in the Halo 
Abundance Matching that TGI 1 use - both types of galaxy are 
assumed to have the same baryon mass fractions. Further, gas 
mass fractions are very much lower in ellipticals than they are in 
spirals (star formation has been completed in the former but is 
ongoing in the latter; e.g. see Rodrfgue z-Puebla et al.|20lT[ fig. 
2), yet the empirical M star -M gas relation used in TGI 1 for both 
types of galaxy (from B aldry et al.|2008| l is for spirals only. Thus 
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their theoretical line strictly applies to neither spirals nor ellip- 
ticals on their own, but rather an amalgam of the twoj^] Given 
the uncertainties introduced by ellipticals described in the para- 
graph above (and in light of the considerable amount of con- 
fusion that already exists surrounding the Tully-Fisher Relation 
itself) it seems sensible to remove ellipticals from the discus- 
sion. In the context of HAM, one should therefore match spiral 
galaxies in the observational survey with simulated halos that 
are believed would harbour spirals in reality. The latter may be 
identified by their environment and merger history, and matching 
of this type has, since the time of TGI 1 's study, been performed 
inll 



Rodrfguez-Puebla et al. ( 2011 hereafter RP 1 1 ) 



1.4. Aim 

In view of the above issues, my approach will be the following. 
I aim to construct a theoretical BTFR, taking account, as com- 
pletely as possible, of all relevent effects that are seen in up-to- 
date numerical simulations of dark matter halos. I will use only 
parameters and relations that are expected a priori from analytic 
and numerical studies of dark halos, and will not include any tun- 
able quantities. Where applicable, these relations will be specif- 
ically for spirals. In addition, I will propagate into the BTFR the 
expected scatter in all relations that I use, investigate the effect 
of adiabatic contraction, and determine the robustness of my re- 
sults to changes in the halo density profile. The result will be the 
expected BTFR in the ACDM paradigm, which may be fairly 
compared to a range of observational data sets. Two such data 
sets will be used to qualitatively illustrate the degree of agree- 
ment between theory and observation. Finally, I will investigate 
the extent to which the prediction would be modified by mod- 
elling two selection effects in the observational data. 

The structure of this paper is as follows. In Sect. [2] I describe 
the technical details of my calculation of the predicted BTFR. 
Section [3] explores the issues involved in selecting appropriate 
observational data for comparison with theory, and describes the 
data sets that will be used. The results will be presented and dis- 
cussed in Sect. [4] Section [5] contains a comparison of my find- 
ings with those of three recent studies in the literature. Finally, 
Sect, [(^summarises my conclusions and suggests fruitful avenues 
for further work. 



2. Method 

The basic semi-analytic methodology employed in this work is 
adapted from M98 (Sect. 2). As the starting point, I take a spi- 
ral galaxy with a particular stellar mass. The M v i r -M stSLr relation 
from the Halo Abundance Matching performed in RP11 (eq. 5) 
is then used to calculate the virial mass of the surrounding dark 
halo. The corresponding gas mass is calculated from the M star - 
M gas relation described by RP1 1 eq. 2, in order to derive an an- 
alytic form for the HAM M vil -M bm relation (RP11, fig. 5). The 
total halo mass is then used to derive the virial radius, virial ve- 
locity and scale-length of the halo (assumed for the moment to 
have an NFW density profile; see M98, eq. 2). Given these quan- 
tities, I find a self-consistent solution for the scale-length (/?</) of 
the spiral galaxy's disk (assumed exponential) that is formed by 
baryon collapse, and the total rotation curve. R c { is found via cal- 
culation of the total energy of the halo and angular momentum 



3 Incorrect gas mass fractions will have a far larger effect on the 
BTFR than the LTFR. Part of the reason why this discrepency is not 
considered to be very important by TGf 1 is the subordinance of the 
BTFR to the LTFR in their work. 



of the disk (an integral over the rotation curve), and comparison 
with a dimensionless spin parameter A (see M98, eq. 9). This 
parameter is used because its distribution is well-constrained by 
numerical simulations, a fact that will become important in the 
discussion of scatter in Sect. 4.2 Before this, A will be set equal 



to its mean value of 0.05. The response of the dark halo to disk 
formation follows the adiabatic contraction prescription laid out 
in |Gnedin et al. (201 1 1. Ra depends on the concentration of the 
halo, which will be derived from M v ; r according to the results 



of |Zhaoet al.| (|2009), fig. 16. 

Starting with an estimate for Rj, the adiabatic contraction 
equations may be solvecj^Jusing the Newton-Raphson method to 
determine the final density profile of the dark halo and thence 
the rotation curve (a sum in quadrature of disk and halo con- 
tributions). The latter may be used to find a new estimate of 
the disk's angular momentum, and hence of Rj. This updated 
scale-length is then adopted and the process iterated until the 
Rd value generated from the rotation curve is identical to the 
value that produces that curve. (My convergence requirement is 
that f r [see M98, eq. 29] change by less than 0.1% in the fi- 
nal iteration; reducing this is found to make a negligible differ- 
ence to the results.) The self-consistent rotation curve is then 
sampled at 1000 radii logarithmically evenly spaced between 
0.0001 x R vir and 0.4 x R vlI , and 100 radii between 0.4 x R vlI 
and R v i r , in order to determine the maximum rotation velocity of 
a model galaxy with this mass. Finally, this procedure is repeated 
for 200 stellar masses logarithmically evenly spaced in the range 
1O 5 M < M star < 1O I2 M , allowing the predicted BTFR (maxi- 
mum rotational velocity vs baryonic mass) to be plotted. 



3. Observational data 

There have been numerous studies measuring the baryonic mass 
and rotational velocity of actual spiral galaxies, all producing 
results that vary to a greater or lesser extent from one to an- 
other. Reasons for this include (but are not limited to) differ- 
ences in the following: 1) Wavelength used for luminosity mea- 
surements; 2) Mass-to-light ratio used for converting luminosity 
to stellar mass; 3) Method for accounting for mass of atomic, 
molecular and ionised gas; 4) Method for measuring rotational 
velocity (e.g. using line width or a resolved rotation curve); 5) 
Radius at which the rotational velocity is measured (for rotation 
curve studies); 6) Type of galaxy survey; 7) Criteria for decid- 
ing exactly which spirals should be included. The precise ef- 
fects of variations in these factors are poorly-understood, but are 
certain to produce non-trivial systematic shifts in the resulting 
BTFRs (e.g. see Fo reman & Scott|2012| p. 3 and refs. therein). 
The nature of the observational data used for comparison with 
a BTFR prediction is therefore crucial for determining whether 
or not the prediction appears successful. This makes imperative 
the preliminary establishment of rigorous criteria for selecting 
observational data sets which have maximum compatibility with 
a given theoretical model. 

The first step in this endeavour is to investigate the approx- 
imations and assumptions that go into the theory. The Halo 
Abundance Matching used to derive the M v j r -Mt, al relation con- 
tains both a prescription for calculating stellar mass from lumi- 
nosity, and gas mass from stellar mass. Thus the resulting BTFR 



Gnedin et al. 



1 2011 1 demonstrates that no single set of parameter 



values reliably specifies the effect of adiabatic contraction in all cases, 
and recommends considering a range. All plausible values are found to 
produce almost identical TFRs, so for simplicity I will use A = 1.6, w 
= 0.8 (see Gnedin et al. 2011 eq. 4) throughout this work. 
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prediction should be compared with observational data employ- 
ing the same mass-to-light ratios and techniques for gas mass 
measurement. In addition, I use the maximum rotational velocity 
to describe my model galaxies; observed galaxies should use a 
similar rotation measure. By using only observations with identi- 
cal assumptions to a theoretical prediction, unknown systematic 
effects associated with variation in the parameters listed above 
will be eliminated. As long as the observations are consistent 
among each other, cuts may be applied to the model galaxies to 
enhance compatibility with the data. In principle (if the observa- 
tional uncertainties were sufficiently well-known) a hypothesis 
test such as x 2 could then provide a quantitative measure of the 
extent of agreement between theory and observation. 

From this perspective, it may be counterproductive to bin 
together the results of many different galaxy surveys without re- 
gard of the different assumptions and approximations on which 
they are based (as done for example in TG11, fig. 11). The ad- 
vantage of comparison with a large quantity of observational 
data is offset by a reduction in quality in the sense that the ob- 
servational details of the surveys are lost. Further undesirable 
features of this method of data presentation are the following: 

1 . The error bars show the scatter of the data points in each bin 
around their mean. Information concerning the uncertainties 
on the individual measurements is entirely lost, making it 
impossible to judge whether or not the points could actually 
be consistent with the theoretical line. 

2. It is not known how many data points are in each bin. Some 
bins may contain many more than the others, and should 
therefore be weighted more heavily when comparing to the 
theoretical line. This cannot be taken into account in a visual 
assessment. 

3. The mass and velocity range of the data points in each bin 
is unknown - the standard deviation is the only measure of 
dispersion that is retained. 

4. Projecting each data point onto one of a discrete set of rota- 
tional velocity values causes loss of accuracy. 

I will use two observational data sets for comparison with the 
BTFR prediction generated by the methodology of Sect. [2] I do 
not intend this to constitute a definitive test of theory, but rather a 
rough qualitative assessment of my results: a comprehensive lit- 
erature review to find all data sets consistent with the principles 
above is beyond the scope of this paper. Indeed, my intention is 
to produce a prediction that includes a minimum of theoretical 
assumptions and may therefore be compared to a range of data 
sets as others see fit (and indeed as BTFR observations improve 
and proliferate). In the next two sub-sections I will describe the 
data that I use, and my reasons for selecting them. 

3.1. Gas-rich spirals 



As discussed in Sect. |1.3.1| a major (if not the dominant) source 
of uncertainty in BTFR measurements is the mass-to-light ra- 
tio used to determine M star . Thus it makes sense to prioritise 
gas- as opposed to star-dominated spirals, for which the rela- 
tive contribution of the uncertainty in M stal to the error budget 
is low. MG12 presents a compilation of 47 spirals which have 
M gas > M star and which moreover satisfy several quality cri- 
teria including a resolved rotation curve and consistent optical 
and HI inclinations. McGaugh uses the flat part of the rotation 
curve as the rotational velocity measure - this can be approx- 
imately accounted for in my modelling by discarding galaxies 
with non-flat rotation curves (see Sect. 4.3 1. A further discrep- 
ancy between McGaugh's approximations and those of Sect.|2]is 



that RP1 l's Halo Abundance Matching uses mass-to-light ratios 
derived from the Yang et al. (2009) stellar mass function, whilst 
McGaugh assumes a Portinari population synthesis model and 
Kroupa Initial Mass Function. However, the gas-dominated na- 
ture of the MG12 galaxies renders them relatively insensitive to 
variations in the mass-to-light ratio. 

3.2. The GHASP survey 

No gas-dominated spirals exist with Mb ar <: 1O 1O ' 5 M . To extend 
the BTFR to higher-mass (e.g. to better constrain its slope), the 
data from Torres- Flores et aL] ( |201 1) will be used. This is de- 
rived from a homogeneous galaxy survey (GHASP) undertaken 
using a scanning Fabry-Perot interferometer in France. Besides 
homogeneity, the GHASP survey offers the advantages of a com- 
pletely resolved 2D velocity field (enabling rotation curves to be 
determined without uncertainties concerning position angle or 
inclination), observation of close galaxies (minimising distance 
uncertainties), and the adoption of V max as the rotational veloc- 
ity measure (in accordance with my theoretical model). Mass-to- 
light ratios are obtained using a Bell population synthesis model, 
which is shown in RP1 l's fig. 1 to produce a similar galaxy stel- 



lar mass function to the Yang et al. (2009) results used in the 
Halo Abundance Matching. 



4. Results and discussion 

4.1. The baryonic Tully-Fisher Relation predicted by ACDM 

In this subsection, scatter in the BTFR prediction is entirely sup- 
pressed. Fig. [T] compares three theoretical BTFRs with the ob- 
servational data described in Sect. [5] (recQ and blue points)]^] In 
each case, an NFW halo density profile is assumed. The black 
line uses the M v j r -Mb ar relation for spirals only (RP11 table 1, 
column 3) and includes the effect of adiabatic contraction, and 
is therefore the primary result of this subsection. This predic- 
tion is contrasted with two others, one for which adiabatic con- 
traction is switched off (i.e. halo unaffected by disk formation; 
cyan) and another which uses the M v j r -Mb ar relation for ellipti- 
cals (magenta). This takes M v ; r to be calculated from M star using 
the results in column of 4 RPll's table 1, and M gas from the 
red line in RP11, fig. 2, and is displayed to illustrate the extent 
to which contamination of the HAM M v i r -M\, ar relation with el- 
liptical galaxies may be expected to distort the BTFR (see also 
Sect. 5.1 1. Some observations from this plot: 



1 . Each line is almost a power- law up to M bar ~ lO la5 M , be- 
yond which it becomes much shallower. This is clearly due 
to the nature of the M^-M^ relation, in which a similar 
break occurs (RP11, fig. 5, top right panel). At high virial 
masses, M^ m rises very slowly with M v ;,. Thus the increase 
in rotational velocity is large relative to that in baryonic mass 
and the BTFR flattens out. 

2. The curvature becomes very obvious around M\, m w 10 n M G . 
However, this is approximately where the observational data 
ends. The data does not, therefore, provide compelling evi- 
dence for either a continuation of the power-law from lower 
masses (as expected for example in MOND) or for the de- 
crease in slope predicted here. Observations of higher-mass 



5 One of the galaxies from MG12 (DDO 210) was removed due to its 
extremely low baryonic mass. 

6 The green data points are from TGI 1. A comparison with these will 
be presented in Sect. 5.1 
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Fig. 1. The black (magenta) line uses the M v j r -M star and M star - 
M gas relations for spirals (ellipticals) - see RP11, eq. 5 and fig. 
2. The cyan line is the analogue of the black line but with AC 
switched off. Red data points are from McGaugh ( 2012)), blue 
from|Torres-Flores et al. ( 2011| l and green from Trujillo-Gomez 
et al.| ( 201 l| l. The M vlI -M\, m relation must be extrapolated below 
the dotted horizontal line. 



spirals would be needed to decide between these two possi- 
bilities. 

3. There is quite good agreement with the entire GHASP sam- 
ple. However, the predictions lie systematically below the 
data for the gas-dominated spirals. In this respect, the black 
line gives a better fit than the magenta line, suggesting that 
a failure to discriminate between spirals and ellipticals will 
reduce agreement with the data at low baryonic mass (as 
might be expected from the fact that the observations are 
specifically of spirals)]^] It is crucial to note, however, that 
the Myir-Mbai relation from RP11 must be extrapolated for 
Mbm- ^ 1O 8 9 M (see RP11, fig. 5), increasing the uncertainty 
in the theoretical prediction in this region. This threshold is 
shown by the horizontal dotted line. 

4. The effect of AC is to slightly reduce Mb ar for spirals with 
10 2 km s < V rot < 10 2 ' 5 km s . This may marginally im- 
prove agreement with the observational data. 

Fig. [2] shows the results when the density profile of the 
model halos is switched from NFW to E inasto (|Gao et al-12 008 ), 
Moore ( |Moore et al.|1999] >, or Burkert ( |Burkert|1995| l. To con- 
struct these lines, eqs. 19, 21 and 23 of M98 were modified ac- 
cording to the change in p(r). For the Moore and Burkert pro- 
files, an identical M™ -concentration relation to the one used for 
NFW was adoptedF] Concentrations for the Einasto profile were 
derived from M v ; r using eq. 6 of |Gao et al.| ( |2008] ), a study work- 
ing explicitly with this density profile. The Einasto shape pa- 
rameter a was calculated using |Gao et aL] ( |2008| l, fig. 2. The 
curvature in the BTFR prediction is largest using the Einasto 
profile and smallest using Burkert, but, in general, changing the 
density profile has relatively little effect. Whilst the Einasto re- 
sult fits the data somewhat better at the highest and lowest bary- 



Fig.2. Data points as in Fig.[l] See text for further details. 



onic masses, the NFW line appears near-optimal in the range 
1O 9 M <M bar < 10 n M o . 



4.2. Intrinsic scatter 

Three sources of scatter in the theoretical BTFR are variability 
in halo spin, concentration and baryon mass fraction. The spin 
parameter A is determined by cosmological torques and inter- 
galactic tidal interactions ( |Eisenstein & Loeb|1996| l, and there- 
fore varies with environment. Numerical simulations find A fol- 
lowing a log-normal probability distribution with mean 0.05 and 
standard deviation 0.5 (Mo et al. 1998). Galaxies of the same 
mass may have different concentrations (differently shaped den- 
sity profiles) due to differences in their merger history and hence 
epoch of virialisation. Although no uncertainty is quoted in the 
source of the M^-concentration relation used here (Zha o et al. 
2009), other studies have indicated that the concentration fol- 
lows a log-normal distribution with standard deviation roughly 
independent of M v ; r and equal to 0.18 (Jing 2000; Bullock et al. 
|2001| . Scatter in the HAM M v ; r -M bar relation comes from a vari- 
ety of sources (RP11, Sect. 3.1.1; see also Behroozi et al. 2010). 



RP1 1 quote the scatter in log M bar as 0.23 dex at all virial masses. 

To model the effect that these sources of scatter have on the 
intrinsic scatter of the BTFR, I make the following modifica- 
tions to the methods used to generate the black line of Fig. [T[ 
1) The number of M star values is increased from 200 to 500. 2) 
For each M star , 500 different values of A, M bar and concentration 
are randomly drawn from their respective probability distribu- 
tions. The mean of the log-normal M bar and concentration dis- 



5.1 



7 This issue will be discussed further in Sect. 

8 Concentration depends on halo scale-length, which strictly is only 
defined for the NFW profile. For alternatives, the scale-length will be 
taken as the radius at which the density becomes proportional to r~ 2 . 



tributions are the values of these parameters used in Sect. 4.1 
3) A 150x150 element grid is constructed spanning the range 
10 6 Mq < M bar < 10 12 M Q , lO 07 km s" 1 < V max < 10 2 7 km s" 1 . 
For each of the 500x500 = 250 000 input parameter sets, M bar 
and V max are calculated and identified with a particular element 
of this grid. The final number of points in each element is then 
outputted and colour-coded for comparison with the observa- 
tional data. 

The resulting contour plot constitutes Fig. [3] and shows the 
band in which observational points are predicted to lie, taking 
into account the expected intrinsic scatter. A little under 250 000 
model galaxies lie within the mass and velocity ranges shown in 
this figure, and a grid element is coloured if it contains at least 
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Comparison of data with BTFR including intrinsic scatter 



]_-]_ 5 McGaugh Data 
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Predicted BTFR band with baryon angular momentum halved 
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Fig. 3. Dependence of model galaxy count on region of the M^. dr - Fig. 4. As Fig. [3] but with the baryon angular momentum frac- 
Vmax plane, following the procedure of Sect. 4.2 The contours tions of all model galaxies halved. 



are at 10, 50, 100, 150, 200, 250, 300, 350, and 400. 



10 galaxy points. Thus it should be rare, in the ACDM model 
presented here, for galaxies to lie outside the coloured band. As 
with the black line in Fig. [T] there is good agreement with the 
data at Mb ar ^ 10 9 M Q and poor agreement below (where the 
HAM M v i r -Mbar relation must be extrapolated). The amount of 
intrinsic scatter looks to be in reasonable accord with that in the 
data. 

An additional quantity required for calculation of the maxi- 
mum rotation velocity of each model galaxy is the fraction of the 
galaxy's total angular momentum that belongs to the baryons 
(jd', see M98, eq. 8). Thus far, jd has been set equal to M bar / 
M„j r , the fiducial assumption of M98. However, this paper also 
cites evidence from numerical simulations (in Sect. 2.2) that this 
may overestimate jd- To test the robustness of my results to a 
decrease in baryon angular momentum fraction, I plot in Fig [4] 
the analogue of Fig. [3] but with all jd values halved. The effect 
is seen to be a slight downward shift of the predicted band at 
A^bar ^ 1O 9 ' 5 M (i.e. increased curvature) in addition to a small 
increase in scatter in this region. This decreases agreement with 
the high-mass data points, but to a sufficiently marginal extent 
for plausible changes in jd to be considered inconsequential. 

4.3. Selection effects 

In Sect. |4.2| all allowed halo spins, concentrations and baryon 
mass fractions were used to construct model galaxies, and con- 
tributed to Fig. [3] However, it is unlikely that all such parameter 
values would yield spiral galaxies whose properties astronomers 
would measure and plot on the BTFR. Reasons for this fall in 
one of two categories: 1) Halos with extreme properties may be 
unstable and hence never form disk galaxies; 2) Astronomers 
typically require galaxies to fulfill certain selection criteria to 
be included in the TFR, with the intention of minimizing sys- 
tematic uncertainties and ensuring some degree of consistency 
within the sample. A proper comparison of a theoretical BTFR 
with observational data should limit the input parameter space to 
the regions which generate stable spiral galaxies with properties 
passing all selection requirements of the data (see also Sect. [3}. 
In this section, I explore two such selection effects, one from 
each of the categories listed above. 



It is theoretically expected that low A will make a spiral 
galaxy prone to the development of a bar instability, potentially 
causing either the total disruption of the disk or its transforma- 
tion into an irregular galaxy which would not be plotted on the 
TFR (see e.g. Christodoulou et al. 1995] >. Removing the unstable 
model galaxies from the contour plot may therefore boost com- 
patibility with the observational data. A simple prescription for 
the stability threshold is given by M98 (Sect. 3.2) in terms of 
the mass and scale-length of the disk, and maximum rotational 
velocity of the halo. Since these quantites are calculated for each 
model galaxy as part of the BTFR calculation, removing unsta- 
ble galaxies in this way is straightforward. 

Of the ~250000 model galaxies, only 5 169 are rejected us- 
ing the stability threshold e mfirll = 1, 10475 for e mj „.„ = 1.1 and 
17 512 adopting the the upper limit deemed plausible by M98, 
e m ,crit = 1-2 (a larger value causes rejection of halos with larger 
spins; see M98 eqs. 35 and 37). Thus we expect that the effect on 
the BTFR of rejecting galaxies unstable by this criterion is small, 
indicating that most of the model galaxies are in fact stable. This 
is illustrated in Fig. [5] a replica of Fig. |3]but with galaxies un- 
stable to bar formation at the e mfirlt = 1.2 level removed. The 
majority of unstable galaxies have baryonic masses in the range 
10 10 - 10 105 M o , and the theoretical intrinsic scatter in this re- 
gion is minorly reduced by excluding them. At lower-masses, 
where disagreement with the data is most pronounced, there is 
virtually no effect. 

An important observational issue is the point on a galaxy's 
rotation curve at which velocity is measured. MG12 adopts Vfl at 
(due to earlier work suggesting that this minimizes the scatter 
in the TFR), and defines a rotation curve as flat if the difference 
between the velocity at 3 disk scale-lengths and at the last mea- 
sured point of the rotation curve is less than 15% (Stark et al. 
2009 ). This requirement is applied to many galaxies in the gas- 
rich data set, which are therefore not directly comparable to the 
data points for model galaxies with parameter values such that 
their rotation curves do not plateau. Here I apply a somewhat 
modified requirement to remove such galaxies, which should 
nevertheless produce similar results: the difference between the 
maximum velocity and that at 3 disk scale lengths must be less 
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Predicted BTFR band excluding unstable galaxies 




1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 



loglO ( Rotational Velocity / km s 1 ) 

Fig. 5. As Fig. [3] but excluding model galaxies with 
V max / iGM bai /R d < 1.2. See text and M98 (Sect. 3.2) for fur- 
ther details. 

than 15%|^] Only 200 model galaxies are excluded for having 
non-flat rotation curves, resulting in a negligible change to the 
BTFR band. Since rotation curves are typically observed to be 
flat in the outer regions of spiral galaxies, this gives some confi- 
dence that the galaxies produced by the methods and relations 
of Sect. |2| ar e realistic. These results also support the claim 
of |McGaugh| ( [20T2) and |Stark et aLlp009) that V flat and V max 
are typically very similar. 

5. Comparison with the literature 

In this section I relate my results to those of TG11, D07 and 
MG12. 

5. 1 . Trujillo-Gomez et al. (2011) 

TG11 use the Bolshoi N-body simulation in association with 
baryon mass fractions from HAM to give a prediction for the 
BTFR which exhibits considerable curvature (their fig. 11). This 
presumably derives from the precise M v j r -M bal - relation that is 
adopted (although this is not actually visible in TG11). The 
fact that no such curvature is evident in the data is pointed out 
in MG12, where it is used to quickly dismiss TGll's predic- 
tion; yet TG11 claim to find a "good fit" to the observational 
data. The extent to which curvature militates against a predicted 
BTFR depends on the status that the relation has within the 
paradigm under consideration. For TG11, the BTFR is simply 
an observed empirical relation that arises in some complicated 
way via galaxy formation, growth and merger history in addition 
to baryonic processes. The result is that, from this perspective, 
there is no reason for the BTFR to be linear in log-log space, and 
it is a coincidence if it approximately is. As long as the curved 
prediction is more or less consistent with each data point, the 
agreement may be considered adequate. The situation is very 

9 The reason why McGaugh's flatness criterion cannot be applied 
directly is that the "last measured point" is not defined for the model 
galaxies. To improve the theoretical flatness requirement, one would re- 
quire a radius marking the end of the observationally-resolved rotation 
curve. This radius would depend on the observational techniques and 
instruments used in the galaxy survey under consideration. 



different, however, for McGaugh, who holds the BTFR to be a 
relation arising directly from a fundamental force law (MOND). 
From this persepective, the perfect linearity of the BTFR is one 
its defining characteristics, and hence inability to reproduce this 
is a major shortcoming. Furthermore, it can be seen in Fig. [T] 
(or by comparison of MG12's fig. 5 with TGll's fig. 11) that 
although a power-law fit to both MG12's and TGll's data sets 
would have a best-fit slope of around 4, the intercept in the TGI 1 
data is significantly lower. This increases agreement with their 
prediction. 

Aside from the issue of whether or not the linearity of the 
BTFR is one of its fundamental properties, it is clear from fig. 
1 1 in TGI 1 that the data points do lie above the predicted line at 
both low and high masses: the curvature does appear to be too 
great. This may be partially explained with the help of Fig. [T 
in which the magenta line shows the analogue of the theoretica 
BTFR for elliptical galaxies. Contamination of TGll's HAM 
M v i r -M star relation with ellipticals would be expected to give a 
result between the red and magenta lines, which might exhibit 
considerable curvature p"| One might also expect to be able to 
explain the behaviour at the high-mass end of fig. 1 1 in this way. 
TGll's prediction lies below the data points describing spirals 
but is in agreement with those describing ellipticals. However, 
my Fig. [T] shows that exchanging the expected spiral baryon 
mass fractions with those for ellipticals has very little effect on 
the predicted BTFR at high mass (if anything, the relation for el- 
lipticals lies above that for spirals). A possible explanation is that 
not only does the M v j r -M bal - relation depend on galaxy morphol- 
ogy, but other relations important for calculation of the BTFR 
do too (e.g. the M vu -concentration relation). This reinforces the 
point made in Sect. |1.3.2| that the differences between the fun- 
damental parameters of spirals and ellipticals are not sufficiently 
well-known to permit a fair comparison between them in bary- 
onic mass-rotational velocity spacef^j 

A further difference between the Abundance Matching of 
TGI 1 and RP1 1 is that the former uses only the galaxies in their 
TFR sample, whilst the latter matches the theoretical halo mass 
function to the entirety of the Sloan Digital Sky Survey. Thus the 
Mvir-M s tar relation used by TGI 1 is built from around 1 000 data 
points, whilst that described in RP11 contains of order 500000 
and should therefore enjoy greater precision. Underneath this 
difference runs the question of whether the BTFR is a funda- 
mental relation applicable to spirals of any mass, composition 
and location (a view championed by McGaugh), or whether it is 
a product of complex evolution histories in a ACDM universe. 
In the latter scenario, the BTFR may well be different for dif- 
ferent sets of galaxies - which may have been subject to signif- 
icantly different conditions during their evolution - and hence 
it would be reasonable to limit the M v i,--M bal - relation to those 
galaxies used in the TFR. However, if all types of spiral lie on 
the same BTFR, it is unsatisfactory to "explain" the BTFR for 
only a small subset of these. A universal M v i r -M bar relation (as 
approximated for example by the results of RP11) is necessary 
to give the theoretical BTFR universal applicability. 



10 However, the fact that spirals are far more abundant than ellipticals 
at low masses suggests that the HAM Mvir-Mbar relation for all galaxies 
- as used by TGI 1 - should be much closer to the spiral line than the 
elliptical line. 

11 Of course, the theoretical difference between spirals and ellipticals 
in this space is known only to the extent that the halos hosting spirals 
and ellipticals may be accurately distinguished (as described in RP11, 
Sect. 2.2. 1). Methods for achieving this are still in their infancy and may 
be significantly refined in the future. 
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5.2. Dutton et al. (2007) 

D07 use a semi-analytic approach, also based on M98, to 
investigate the LTFR in the ACDM paradigm, using simple 
parametrizations of the required relations which include tunable 
parameters. These are selected to fit the size-lumunosity (RL) 
relation in addition to the LTFR. D07 find that the slope and 
intercept of the M v j,-Mb ar relation can be tuned to fit the RL re- 
lation without significantly affecting the TFR. While it may be 
true that baryon mass fractions are of little importance to the lu- 
minosity TFR (the primary concern of D07), they are clearly of 
great importance to the baryonic TFR. For example, the BTFR 
described by the black line in Fig.[T]is significantly steeper than 
early predictions for the BTFR adopting fixed baryon fractions 
(eg. Mo et al.[T 998 1 Bullock et al. 2001) primarily because of the 
steep M v i r -Mbar relation predicted by HAM. Thus we see that 
successful prediction of the BTFR does not necessarily follow 
from successful prediction of the LTFR (cf. Sect. 



ing of the way in which the appropriate number of baryons is 
ejected during galaxy formation (and where they went) remains 
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In a similar vein, D07 claim that scatter in baryon mass frac- 
tion does not significantly affect the scatter in the LTFR, and, in 
view of the over-budgeted scatter in the RL relation, they set this 
to zero. Not only does this lack physical motivation (if the low 
observed baryon mass fractions are due to expulsion of baryons 
by astrophysical phenomena, one would clearly expect statistical 
fluctuations between galaxies), but it again neglects the fact that 
scatter in baryonic mass has a large effect on the BTFR. Indeed, 
this contributes a significant fraction of the scatter in Fig. [3] The 
results in D07 are specific to the luminosity TFR, and many of 
the conclusions are unlikely to hold if the BTFR were considered 
also. 

Finally, D07 argue that adiabatic contraction must be re- 
placed by expansion in order to successfully reproduce the LTFR 
and galaxy luminosity function with a realistic Initial Mass 
Function. Although I make no claim to "successfully reproduce" 
the BTFR data with my model, I do notice from Fig. [T]that the 
effect of adiabatic contraction is very small and does, if anything, 
increase agreement with the data. D07's result may be specific 
to the parametrizations and parameter values produced by their 
particular optimisation procedure. 

5.3. McGaugh (2012) 

McGaugh argues that the failure of ACDM models to repro- 
duce the low intrinsic scatter and high slope (~4) of a power-law 
fit to the BTFR without excessive fine tuning constitutes evi- 
dence against the paradigm. Given the existence of a Modified 
Gravity Model (MOND) which naturally explains these features, 
McGaugh further claims that the BTFR militates against the ex- 
istence of dark matter itself. In this subsection I discuss the way 
in which my results impinge on two aspects of this argument. 

In MG12, McGaugh claims that natural ACDM predictions 
have a slope of approximately 3. This was indeed the case for 
the "naive" derivation of Sect. 1 1 . 1 1 in addition to several subse- 
quent studies that introduced a dependence of halo concentra- 
tion on M v ; r (e.g. |Bullock et al.|[200T| ). However, the slope is 
significantly steepened by replacing the constant baryon mass 
fractions used in these works by more "realistic" values, such as 
those derived from Halo Abundance Matching. For example, the 
low-mass end of the black line in Fig.[T]has a slope in excess of 
5. On the one hand, the HAM M v h--Mb ar relation could be said 
to be a "natural" ACDM prediction in the sense that the halo 
mass function used in the Abundance Matching is determined 
from N-body simulations of dark halos in a ACDM cosmolog- 
ical background. On the other hand, a theoretical understand- 
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is a more fundamental problem for ACDM than failure to accu- 
rately reproduce the observed BTFR. Indeed, it is possible that 
the remaining discrepency between the observed and predicted 
BTFRs in Figs. [T] and [3] is entirely attributable to errors in the 
adopted M^-M\,. M relation. It is clear that this relation makes 
a huge difference to the BTFR (contrast my results with those 
of TGI 1 and |Bullock et al.pXH) , yet the method of HAM is 
only just now being fully developed, and may be subject to ma- 
jor refinements in the future. We also see from Fig. [3] that the 
BTFR prediction is significantly different to the data only in the 
region in which the HAM M v j r -Mb ar relation must be extrapo- 
lated. Could this be due to the form of this relation changing at 
low halo masses? 

McGaugh claims that the scatter of gas-rich galaxies around 
a power-law BTFR may be fully accounted for by the uncer- 
tainties on the data points. Even if this were soj^jit would not 
preclude the possibility of a non-zero intrinsic scatter in the un- 
derlying theoretical relation. Many potential BTFR predictions, 
with a range of non-zero intrinsic scatters, may be "compatible" 
with the data in the sense that they cannot be rejected at the 95% 
confidence level. McGaugh's argument succeeds in demonstrat- 
ing that the data is consistent with MOND (which no doubt is 
worthy of consideration for providing a simple explanation of 
this and many other aspects of galaxy phenomenology), but not 
that it is inconsistent with ACDM. Indeed, the magnitude of the 
theoretical scatter in Fig. [3] does not appear wildly discrepant 
with that in the data. 



6. Summary and suggestions for further work 

This paper illustrates methodically the way in which the "naive" 
slope-3 ACDM prediction for the BTFR can be improved by 
taking proper account of concentration, baryon mass fraction, 
and baryon disk collapse followed by adiabatic halo contraction. 
The slope of the relation is significantly increased by adopt- 
ing baryon mass fractions from spiral-only Halo Abundance 
Matching, improving agreement with observational data at mod- 
erate to high baryonic masses. Predicted rotational velocities are 
systematically higher than those observed at the low-mass end, 
although plotting the theoretical BTFR in this region requires ex- 
trapolation of the Myir-Mbar relation obtained by the Abundance 
Matching. In addition, curvature (in log-log space) is introduced 
into the prediction and becomes significant for Mb ar > 10 n M G . 
Although curvature is not evident in the data, spirals with such a 
high baryonic mass are rare, and the observational studies used 
here do not preclude a reduction in the slope of the BTFR at 
very high masses. The effect of adiabatic contraction is small, 
especially at low and very high baryonic masses. The BTFR pre- 
diction obtained here is reasonably robust against changes to the 
density profile used to describe the dark matter halos, and the 
angular momentum fraction posessed by the baryons. 

Accounting for the expected variability of the halo spin pa- 
rameter, concentration, and baryon mass fraction, I find the the- 
oretical BTFR to have an intrinsic scatter that is not clearly dis- 
crepant with that in the data, although may be somewhat larger 



12 Significant controversy surrounds this issue. For example, Foreman 
|& Scott| J20121 find McGaugh's gas rich data to be inconsistent with 
any zero-intrinsic scatter relation at the 95% confidence level, and many 
other observational studies of the BTFR find at least some evidence for 
a non-zero intrinsic scatter. 
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than necessary to successfully account for the observations. Two 
simple selection requirements imposed on the model galaxies - 
a flat rotation curve and stability against bar formation - were 
found to have minimal effect on the predicted BTFR. This sug- 
gests that the spirals formed in the ACDM model presented here 
have rotation curves that plateau (as typically observed) and are 
predominantly stable. 

There are two main ways in which this study could be use- 
fully supplemented. 

1. Only a subset of the effects causing intrinsic scatter in the 
TFR were considered here. Further sources of scatter in- 
clude the possibility that dark halos have yet to virialise or 
are not well-fitted by an NFW profile (Jing 2000), varia- 
tions in galaxies' mass aggregati on histories (|Eisenstein & 
Loeb 1996}, and halo triaxiality (Kazantzidis et al. 2010] >. 
However, there exist also unconsidered selection factors 
which will reduce the scatter, such as the possibility that disk 
formation in halos with high baryon mass fraction is likely 
to create elliptical or SO galaxies to which the TFR does not 
apply (M ayer & Moore| 2004[ l. The true theoretical intrin- 
sic scatter can only be determined once all of these effects 
have been included in the style of Sect. |4.3| Proper imple- 
mentation of these selection criteria is contingent on a solid 
theoretical understanding of galaxy formation. 

2. The comparison with observational data presented in Sect. [4] 
is far from optimal. As discussed in Sect. [3] the ideal 
would be to amalgamate all data sets that employ compati- 
ble mass-to-light ratios, prescriptions for gas mass measure- 
ment, methods of rotational velocity measurement and selec- 
tion criteria. By implementing similar constraints in the the- 
oretical model, a comparison could be made between theory 
and observation that would be free from systematic errors. A 
simple hypothesis test could then be used to give a quanti- 
tative measure of the agreement between the BTFR and the 
ACDM paradigm as we currently understand it. 

The results of this study will be modified by advances in 
the techniques of Halo Abundance Matching and halo spin and 
concentration estimation from N-body simulations. The modular 
approach will make such updates easy to implement, and their 
effects on the BTFR readily visible. 
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